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A  geometrical  theory  is  devised  for  the  description  and  calculation  of  surface  waves,  such  as  Rayleigh 
waves,  on  boundaries  or  interfaces  of  elastic  solids.  It  applies  to  curved  surfaces  and  to  inhomogeneous 
media.  The  theory,  which  is  an  extension  of  geometrical  optics,  involves  complex  rays  that  travel  from  the 
source  to  the  surface,  then  along  the  surface,  and  finally  from  the  surface  to  points  in  the  solid.  It  also  in- 
cludes phases  and  amplitudes  associated  with  each  point  on  each  ray.  Geometrical  formulas  are  derived  for 
the  determination  of  these  phases  and  amplitudes.  They  involve  certain  excitation  and  radiation  coefficients, 
which  are  also  determined.  The  total  displacement  at  a  point  is  the  sum  of  the  displacements  on  all  the  rays 
through  the  point,  each  of  which  is  constructed  from  the  corresponding  phase  and  amphtude.  The  theory 
applies  to  periodic  waves  of  high  frequency  and  short  wavelength,  as  well  as  to  the  rapidly  varying  portions 
of  any  waveform.  It  is  a  generalization  of  the  authors'  earlier  geometrical  theorv  of  scalar  surface  waves 
[J.  Appl.  Phys.  31,  1039  (I960)]. 


INTRODUCTION 

A  SURFACE  wave  is  a  wave  that  travels  along  a 
surface  near  which  the  wavemotion  is  confined, 
e.g.,  a  Rayleigh  wave.  Because  of  the  difficulty  of 
solving  exactly  elasticity  problems  involving  such 
waves,  we  have  devised  a  geometrical  method  for  their 
approximate  solution.  It  is  an  extension  of  our  theory  of 
scalar  surface  waves'  and  is  based  upon  geometrical 
optics.  Therefore,  it  also  provides  a  geometrical  picture 
of  how  surface  waves  propagate.  It  is  vahd  for  periodic 
waves  of  high  frequency  and  short  wavelength  and  for 
the  rapidly  varying  portions  of  any  waveform.  How- 
ever, experience  with  other  problems  of  wave  propaga- 
tion shows  that  it  is  accurate  even  for  wavelengths  as 
large  as  other  dimensions  in  the  problem. 

Geometrical  optics  determines  the  rays  along  which 
waves  travel.  It  also  yields  the  phase  or  arrival  time  at 
each  point,  and  can  be  used  to  find  the  amplitude  and 
direction  of  the  wavemotion.  This  leads  to  an  expression 
for  an  elastic  wavemotion  that  is  the  first  term  in  an 
infinite  series.  This  has  been  shown  by  V.  M.  Babich  and 
A.  S.  Alekseev-  and  by  us,  and  we  have  shown  how  to 


'  J.  B.  Keller  and  F.  C.  Karal,  Jr.,  "Surface  Wave  Excitation 
and  Propagation,"  J.  Appl.  Phvs.  31,  1039-1046  (1960). 

'V.  M.  Babich  and  A.  S.  Alekseev,  "The  Ray  Method  for 
Calculating  the  Intensity  of  Wave  Fronts,"  Bull.  Acad.  Sci.  USSR, 
Geophys.  Ser.  No.  1,  9-15  (1958). 


determine  the  subsequent  terms.'  However,  although 
this  series  improves  the  description  of  the  motion  on 
the  usual  ray,  it  fails  to  describe  dilTracted  and  surface 
waves.  An  additional  series  is  required  for  each  such 
wave.  We  have  shown  how  to  construct  the  leading  term 
in  each  series  by  an  extension  of  geometrical  optics, 
called  the  "geometrical  theory  of  diffraction."''  That 
theory  has  been  applied  successfully  to  various  scalar-, 
acoustic-,  and  electromagnetic-diffraction  problems,  and 
F.  Gilbert''  has  applied  it  to  certain  elastic-diffraction 
problems.  Now,  we  show  how  to  apply  it  to  elastic 
surface  waves. 

In  the  next  section,  we  introduce  various  new  rays. 
In  Sec.  II,  we  determine  the  phase  and  amplitude  of  the 
wavemotion  on  each  ray  in  terms  of  certain  excitation 
and  radiation  coefficients.  These  coefficients  are  pro- 
portional to  each  other,  as  we  show  in  Sec.  Ill  by  means 
of  the  reciprocity  principle.  In  Sec.  I\',  we  apply  the 
method  to  determine  the  Rayleigh  wave  produced  on 


3  F.  C.  Karal,  Jr.,  and  J.  B.  Keller,  "Elastic  Wave  Propagation 
in  Homogeneous  and  Inhomogeneous  Media,"  J.  Acoust.  Soc. 
Am.  31,  694-705  (1959). 

■•J.  IB.  Keller,  "A  Geometrical  Theory  of  Diffraction,"  in 
Proceedings  of  Symposia  in  Applied  Mathematics,  Volume  8: 
Calculus  of  Variations  and  Its  Applications,  edited  by  L.  M. 
Graves  (McGraw-Hill  Book  Co.,  Inc.,  New  York,  1958). 

'  F.  Gilbert,  "Scattering  of  Impulsive  Elastic  Waves  by  a 
Smooth  Convex  Cylinder,"  J.  Acoust.  Soc.  Am.  32,  841-857 
(1960). 
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the  plane  surface  on  an  elastic  half-space  by  a  compres- 
sional  line  source.  By  comparing  the  result  with  the 
known  short-wavelength  asymptotic  form  of  the  exact 
solution  of  the  problem,  we  determine  the  excitation 
coefficient.  This  coefficient  can  then  be  used  for  any 
wave  incident  on  anj'  surface.  All  the  preceding  results 
apply  to  time-harmonic  waves.  In  Sec.  \',  we  transform 
them  to  apply  to  transient  waves. 

The  theorv',  which  can  be  applied  easily  to  various 
problems,  yields  the  first  term  in  a  series.  Further  terms 
can  be  derived  by  methods  that  have  been  used  in 
similar  scalar  and  electromagnetic  problems. 

1.  COMPLEX  RAYS 

With  each  point  x  of  an  isotropic  elastic  medium,  there 
are  associated  two  propagation  velocities  VdiO  and 
Vs{x),  called  the  compressional  and  shear  velocities. 
In  terms  of  each  velocity,  we  can  define  rays  in  space, 
called  compressional  and  shear  rays,  by  means  of 
Fermat's  principle  or  the  differential  equations  that 
follow  from  it.  Each  ray  can  be  represented  parametri- 
cally  as  a  function  x{s)  of  the  arclength  5.  We  also  define 
complex  rays  of  each  type  as  complex  valued  functions 
x(s)  that  satisfy  the  ray  differential  ecjuations  or  that 
are  determined  by  Fermat's  principle  applied  to  complex 
curves.  This  definition  is  possible,  provided  the  velocities 
are  defined  for  complex  values  of  x,  which  is  the  case  if 
they  are  analj'tic  functions.  Then,  the  complex 
rays  are  curves  in  a  three-dimensional  complex  space. 
Of  course,  our  primary  interest  in  them  concerns  the 
points  where  they  cross  the  real  space.  In  unbounded 
media,  the  complex  rays  that  start  from  a  point  in  the 
real  space  and  then  leave  this  space  usually  never  return 
to  it.  But,  as  we  see  here,  in  bounded  media  the  complex 
rays  may  interact  with  the  boundary  and  then  return 
to  the  real  space,  in  analogy  with  the  fact  that  real 
rays  can  return  to  the  vicinity  of  their  source  after 
reflection  from  a  boundary. 

Let  us  now  consider  a  surface  5  along  which  a  surface 
wave  can  propagate  with  the  velocity  Vr(x).  Real 
surface  rays  on  5  are  curves  that  satisfy  the  differential 
equation  for  rays  traveling  on  a  surface  with  the  velocity 
Vr{x)  or  that  are  determined  by  the  Fermat  principle 
for  curves  on  a  surface.  We  define  complex  surface  rays 
by  applying  Fermat's  principle  or  the  ray  differential 
equations  to  complex  curves  on  .S'.  As  before,  this  is 
possible  if  I'flfx)  and  5  are  defined  in  the  complex  x 
space. 

When  a  space  ray  hits  a  surface  5  it  gives  rise  to  two 
reflected  rays — one  shear  and  one  compressional — and 
two  transmitted  rays  (provided  there  is  a  medium  on 
the  "transmitted"  side  of  S).  The  directions  of  these 
rays  are  determined  by  the  appropriate  laws  of  reflec- 
tion and  refraction.  These  laws  are  consequences  of 
Fermat's  principle  applied  to  curves  having  a  point  on 
the  surface.  It  follows  from  them  that  a  real  incident 
ray  produces  two  real  reflected  rays;  but  one  or  both  of 


the  transmitted  rays  may  be  complex.  The  occurrence 
of  complex  transmitted  rays  is  called  total  reflection  and 
the  waves  associated  with  them  are  called  evanescent 
waves. 

If  surface  waves  can  propagate  along  S,  we  assume 
that  the  incident  ray  may  also  produce  a  surface  ray. 
Then,  the  directions  of  the  incident  and  surface  rays 
are  also  related  by  the  law  of  refraction  in  which  the 
surface-ray  velocity  plays  the  role  of  the  velocity  along 
the  transmitted  ray.  This  is  a  consequence  of  applying 
Fermat's  principle  to  curves  having  an  arc  on  5  on 
which  the  velocity  is  Vr(x)  and  an  arc  in  space.  Since 
the  surface  ray  lies  in  S,  it  follows  from  the  law  of 
refraction  that  the  angle  6  between  the  incident  ray  and 
the  normal  to  5  is  given  by 


sme=V{x)/\'it(x). 


(1) 


Here,  ]'(x)  is  either  l'c(x)  or  r,(x),  according  as  the 
incident  ray  is  a  compressional  or  shear  ray.  This  same 
relation  holds  when  a  space  ray  is  shed  by  the  surface 
ray,  which  we  also  assume  to  occur.  In  addition,  the 
incident-ray  direction,  the  surface-ray  direction,  and  the 
surface  normal  must  all  lie  in  the  same  plane. 

Surface  waves  are  usually  slow  waves,  which  means 
that  they  travel  slower  than  space  waves  and,  conse- 
quently, r(x)/l'ft(x)>  1  at  real  points  x.  But  then,  by 
Eq.  (1),  6  cannot  be  real  at  real  points,  so,  in  this  case,  a 
real  ray  cannot  produce  a  real  surface  ray.  This  is  why 
we  have  to  introduce  complex  rays  to  construct  a  ray 
description  of  surface-wave  propagation. 

There  are  certain  surface  waves  that  travel  with  the 
space  velocity  on  the  convex  side  of  the  surface.  We 
have  previously  called  them  diffracted  surface  waves 
and  the  associated  rays  diffracted  surface  rays.  In 
isotropic  elastic  media,  there  are  two  kinds:  diffracted 
compressional  and  shear  surface  rays.  When  the  incident 
ray  and  the  diffracted  surface  ray  are  of  the  same  type, 
then  r(x)=rfi(x)  and  it  follows  from  Eq.  (1)  that 
6  =  -!r/2.  Thus,  an  incident  ray  of  either  kind  tangent  to 
the  surface  can  produce  a  diffracted  surface  ray  of  the 
same  kind,  and  the  latter  can  shed  space  rays  of  the 
same  kind  along  its  tangents.  W'hen  the  space  ray  and 
the  diffracted  surface  ray  are  of  different  types,  then  it 
follows  from  Eq.  (1)  that  the  diffracted  surface  ray 
must  have  the  greater  velocity,  as  has  been  pointed  out 
by  Gilbert.-^  We  do  not  consider  diffracted  surface 
rays  in  the  present  work. 

Let  us  now  consider  two  points  Q  and  P  in  an  elastic 
medium  bounded  by  a  surface  S  along  which  a  surface 
wave  can  propagate.  There  are  direct  and  reflected  rays 
from  Q  to  P.  In  addition,  a  wave  might  travel  from  Q 
to  P  by  going  from  Q  to  some  point  Qi  on  5,  traveling 
along  5  from  Qi  to  another  point  Pi,  and  then  leaving  5 
at  Pi  to  go  to  P  (see  Fig.  1).  We  suppose  that  this  hap- 
pens, and  we  want  to  determine  Qi  and  Pi.  The  paths 
from  Q  to  Qi  and  from  Pi  to  P  are  space  rays,  while  that 
from  Qi  to  Pi  is  a  surface  ray.  Since  the  rays  QQi  and 
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Fig.  I.  The  complex  ray  from  Q  reaches  the  surface  at  the  com- 
plex point  Qi  where  it  produces  a  complex  surface  ray.  This  ray 
travels  along  the  surface  with  the  velocity  Ffl[x(,5)]  and  sheds 
complex  space  rays.  The  ray  that  leaves  the  surface  at  Pi  passes 
through  P.  Since  the  rays  QQi  and  PPi  may  be  either  shear  or 
compressional  rays,  there  are  four  kinds  of  paths  from  Q  to  P  and 
four  pairs  of  points  Qi,  Pi.  Each  pair  Qi,  Pi  is  determined  by  the 
condition  that  Eq.  (1)  must  be  satisfied  at  each  point  with  the 
appropriate  value  of  F[x(s)]. 

PPi  may  be  either  shear  or  compressional  rays,  there 
are  four  kinds  of  paths  from  Q  to  P  and  four  pairs  of 
points  Qi,  Pi.  Each  pair  Qi,  Pi  is  determined  by  the 
condition  that  Eq.  (1)  must  be  satisfied  at  each  point 
with  the  appropriate  value  of  V. 

11.  PHASE  AND  AMPLITUDE 

We  now  show  how  to  use  the  rays  to  determine  the 
elastic  displacement  u(P)  at  the  point  P.  We  first 
postulate  that  there  is  a  displacement  associated  with 
each  ray  through  P,  and  that  the  total  displacement  at 
P  is  the  sum  of  the  displacements  on  all  the  rays  through 
P.  We  now  assume  that  the  motion  is  time-harmonic 
with  angular  frequency  co,  and  we  write  the  displacement 
as  the  real  part  of  n{P)e~''"',  letting  u  be  complex.  Then, 
on  each  ray  we  write  ii(P)  =  A{P)e"f^^\  and  we  call 
<f(P)  the  phase  at  P  and  A(P)  the  amplitude  at  P. 
This  representation  of  ti  is  made  unicjue  by  the  require- 
ment that  (p{P)  be  proportional  to  cc,  while  the  ampli- 
tude A  may  be  complex  but  may  not  contain  any  phase 
factor  proportional  to  o.\  Finally,  we  assume  that  tp{P) 
is  given  by  the  relation 

viP)='piPi)+w(    ds/Vlx{s)-].  (2) 

J  Pi 

The  integral  in  Eq.  (2)  is  with  respect  to  arclength  5 
along  the  ray  from  Pi  to  P  and  V  is  the  velocity  associ- 
ated with  the  ray,  which  is  either  Vc,  Vs,  or  Vr, 
depending  upon  the  type  of  ray. 

Let  us  apply  Eq.  (2)  to  each  of  the  four  paths  from 
Q  to  P  that  contain  arcs  on  5.  We  may  write  all  four 
results  together  by  introducing  the  subscripts  i  and  j, 
each  of  which  can  take  the  two  values  c  and  s,  denoting 
compressional  and  shear  rays.  Then,  <fij(P)  denotes 
the  phase  at  P  on  a  ray  whose  first  part,  from  Q  to  Qi, 
is  of  type  i  and  whose  last  part,  from  Pi  to  P,  is  of  type 
j.  The  values  of  Qi  and  Pi  on  such  a  ray  are  denoted  by 
Qij  and  P,;.  Then,  we  have 

rQii 


+< 


,f    '  ds/VR\:_xis)-]+o^f     ds/Vlxisn     (3) 
J  Qij  J  Pii 


Here,  *?;(())  denotes  the  phase  at  (3  on  a  ray  of  type  i. 

To  determine  the  variation  along  a  ray  of  the 
amphtude  .1  (P)=  |  A(P)| ,  we  apply  the  principle  of 
energj'  conservation  to  a  tube  of  rays.  For  a  space  ray 
of  type  ;,  this  principle  states  that  A-pVidr  is  constant 
along  the  tube  where  p  is  the  density  of  the  inedium  and 
dr  is  the  cross-sectional  area  of  the  tube.  This  follows 
from  Ecis.  (ILS)  and  (139)  of  Ref.  3  for  compressional 
rays  and  from  Ecj.  (n3)  for  shear  rays.  We  assume  that 
the  same  principle  applies  to  a  strip  of  surface  rays  in 
the  form  A-p'V Rd a  =co\\s,t,  where  da  denotes  the  width 
of  the  strip  and  p  is  a  density,  depending  upon  the 
densities  on  the  two  sides  of  the  surface. 

To  find  the  direction  of  A  (P),  we  utilize  the  facts  that 
on  a  compressional  ray  A  is  tangential  and  on  a  shear 
ray  A  is  normal.  Thus,  in  the  compressional  case,  we 
write  Ac  =  ^ck,  where  k  is  a  unit  vector  tangent  to  the 
ray,  while  in  the  shear  case  we  write  k^  =  Asn,  where  n 
is  the  appropriate  unit  vector  normal  to  the  ray.  Then, 
utilizing  the  energy  principle  that  we  have  on  the  com- 
pressional and  shear  rays,  respectively, 


A,(P)  =  .4,(P,) 
A,(P)  =  vl,(Pi 


■driPi)  p{P,)  F,(P0-l5 


IdriP)    piP)    1 


rfr(Pi)p(Pi)F.(Pi)-|* 
LdriP)    p(P)    ]\(P) 


c(.Pi)T 

k(P),     (4) 

'c(P)J 

]n(P).     (5) 


The  variation  of  n(P)  along  the  ray  depends  upon  the 
torsion  of  the  ray'  and  for  a  straight  ray  n  (P)  is  constant. 
The  direction  of  the  amplitude  associated  with  a  sur- 
face ray  is  completely  determined  by  the  direction  of 
the  ray  and  the  properties  of  the  media  at  the  surface 
and,  therefore,  we  need  not  consider  this  direction.  Thus, 
letting  A  denote  the  amplitude  on  a  surface  ray,  we 
write  merely 


AiP)^A(Pi) 


[dc7\ 
da 


da{Pi)p'(Pi)VRiPi)-ii 
mViP)    VuiP)}' 


(6) 


Let  us  now  determine  the  amplitudes  A,j(P)  on  each 
of  the  four  rays  that  go  from  Q  to  P  hy  way  of  the  sur- 
face S.  To  do  so,  we  make  use  of  Eqs.  (4)- (6)  and  cer- 
tain additional  relations  between  the  amplitudes  on  the 
surface  and  space  rays  at  the  point  on  the  surface  where 
they  meet.  Let  AdQcj)  and  AJQ.j)  be  the  amplitudes 
on  the  incident  compressional  and  shear  rays,  respec- 
tively, at  the  points  Qd  and  ()„■.  Then,  we  assume  that 
the  amplitudes  AdQcj)  and  A,(Qsj)  on  the  resulting 
surface  rays  are  proportional  to  the  incident  amplitudes. 
Thus,  we  write. 

Ac{Qcj)  =  Ec{Qci)AdQcj),  (7) 

A  s  (Q,;)  =  E,  (O.y)  A.  (Qd)  ■  N  (Qd) .  (8) 

In  Eqs.  (7)  and  (8),  EdQcj)  and  EdQsj)  are  excitation 
coefficients  that,  we  assume,  depend  only  upon  the 
properties  of  the  media  at  Q^,  respectively.  In  Eq.  (8), 
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N  (()„•)  denotes  the  unit  normal  to  .V  at  Qsj.  We  have 
introduced  it  because  only  the  component  of  displace- 
ment in  the  plane  of  incidence  is  effective  in  exciting  the 
surface  wave,  and  this  is  proportional  to  the  normal 
component. 

We  next  assume  that  the  amplitude  A:j{]',j)  on  a  ray 
leaving  the  surface  is  proportional  to  the  amplitude 
Ai(P,j)  on  the  surface  ray.  We  express  this  proportional- 
ity in  terms  of  the  radiation  coefficients  Rj{P,j),  which 
depend  upon  the  properties  of  the  media  at  P,j.  Thus, 
we  write 


A,i(l\j)  =  Rj{P,i)A.(P.j)- 


(9) 


When  the  final  ray  is  a  shear  ray,  j=s,  we  assume  that 
the  unit  vector  n (/",>),  which  determines  the  direction 
of  the  displacement  at  Pis,  lies  in  the  plane  containing 
the  tangents  to  the  surface  ray  and  the  shear  ray  at 
P„.  It  is,  of  course,  also  normal  to  the  shear  ray. 
Upon  combining  Eqs.  (4)- (9),  we  obtain 


A.iiP)- 


{ 


dr{P,j)piPii)ViiP.j) 


i{P'i)y 

ViiP)  J 


driP)    piP) 

d<y{Qij)p'iQ.i)  Vnm-^^ 

0- 

XE,iQsj)A,iQ.j).     (10) 


rd<y{Qij)p'iQ.i)  VniQ.i) 

XRiiPj 

ldaiP,i)p'{P,dVn{P,i) 


Let  us  also  write  the  phase  <p,jiP),  given  by  Eq.  (3), 
in  the  form 

+co[     ds/Vlxisn      (11) 

J  p.i 

Then,  we  may  write  the  total  displacement  at  P  due  to 
surface  rays  in  the  form 

a(P)  =  A,AP)Kc{P)e'^"^'''+AUP)K.{P)e'^"'-''' 
+Acs(P)ncs(P)e<^"'''^+A,s{P)nss{P)e'^'-'^K   (12) 

Here,  KciP)  and  k,c(-P)  are  unit  tangents  to  the  com- 
pressional  rays  at  P,  while  n„(P)  and  n„(P)  are  the 
appropriate  unit  normals  to  the  shear  rays  at  P. 

Equation  (12),  together  with  (10)  and  (11),  expresses 
the  main  result  of  our  geometrical  theory  of  elastic 
surface  waves.  It  shows  how  to  express  the  elastic  dis- 
placement at  P  due  to  such  waves  as  a  sum  of  terms. 
There  is  one  term  for  each  ray  through  P  that  comes 
from  a  surface  ray  on  5.  In  Eq.  (12),  we  have  written 
four  terms,  since  usually  there  will  be  four  rays  through 
P  if  the  source  of  the  motion  is  a  point.  For  some  geom- 
etries and  for  other  sources,  there  may  be  more  rays 
and,  therefore,  more  terms  must  be  included  in  Eq.  (12). 
Each  term  involves  the  amplitude  and  phase  on  the 
incident  ray  at  the  point  where  it  produces  the  surface 


ray.  It  also  involves  the  surface-ray  velocity  Vr,  the 
density  p',  the  two  excitation  coefficients  Ec  and  E„ 
and  the  two  radiation  coefficients  Re  and  Rs-  These 
quantities  must  be  determined  by  other  considerations 
— for  example,  from  the  solutions  of  special  or  canonical 
problems.  Once  they  are  known,  the  comuptation  of 
u(P)  is  reduced  to  a  geometrical  calculation. 

In  the  next  section,  we  show  that  the  excitation  coef- 
ficients are  expressible  in  terms  of  the  radiation  coef- 
ficients. In  the  subsequent  sections,  we  determine  these 
coefficients,  as  well  as  Vr  and  p',  for  a  free  surface  by 
using  the  solution  of  a  canonical  problem.  The  same 
method  can  also  be  used  for  any  other  type  of  surface. 

m.  RECIPROCITY  PRINCIPLE 

There  exists  a  relation  between  the  excitation  and 
radiation  coefficients,  as  we  now  show  by  utilizing  the 
reciprocity  principle.'^  This  principle  may  be  formulated 
as  follows: 


"x(/',<3>')  =  "„((?,/'.t). 


(13) 


Here,  u.(P,Qy)  denotes  the  x  component  of  displace- 
ment at  P  due  to  a  unit  force  in  the  y  direction  applied 
to  the  medium  at  Q.  Similarly,  iiy{Q,Px)  is  the  y  com- 
ponent of  displacement  at  Q  due  to  a  unit  force  in  the 
.T  direction  at  P.  We  apply  Eq.  (13)  to  that  part  of  the 
displacement  u  that  involves  surface  waves,  which  we 
have  determined  in  the  preceding  section. 

It  is  simplest  to  apply  Eq.  (13)  to  a  two-dimensional 
problem  in  which  the  media,  the  surface,  and  u  are  all 
independent  of  one  coordinate,  say  s.  The  results  will 
then  apply  only  to  such  two-dimensional  configurations. 
However,  we  assume  that  the  excitation  and  radiation 
coefficients  depend  only  upon  properties  of  the  surface 
and  the  media  in  the  plane  of  incidence,  at  the  points 
of  excitation  and  radiation,  respectively.  Then,  the 
two-dimensional  results  concerning  these  coefficients 
will  apply  in  general. 

On  the  basis  of  these  considerations,  we  consider  the 
displacement  u  produced  by  a  unit  force  in  the  y  direc- 
tion at  Q.  (In  three  dimensions,  this  point  force  repre- 
sents a  line  source  through  Q  parallel  to  the  z  axis.) 
Then,  the  phase  <Pi{Q,j)  can  be  e.xpressed  by  Eq.  (2)  in 
the  form 


V-((?.;)  =  "  ds/Vlx{sn 


(14) 


The  amplitudes  AdQci)  and  A,((),>)  can  be  expressed  in 
the  form  Eqs.  (4)  and  (5).  However,  since  all  the  rays 
lie  in  planes  2  =  const,  the  areas  dr  in  these  and  the  subse- 
quent equations  may  be  taken  to  be  lengths  in  the  xy 
plane  multiplied  by  a  unit  length  in  thes  direction.  For 
the  same  reason,   the  lengths  da  may  be   taken   to 


'  Lord  Ravieigh,  The  Theory  oj  Sound  (Macmillan  and  Co.  Ltd., 
London,  1894;  reprinted,  Dover  Publications,  Inc.,  New  York, 
1945),  Vol.  l,p.  153-155. 
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be    constant    lengths    in    the    z    direction,    so    that 

In  applying  Eqs.  (4)  and  (5),  we  replace  P  by  Q 
and  choose  Pi  to  be  a  point  near  Q.  But  then, 
dT{Pi)  =  S(i6(Q),  where  s  is  the  distance  from  Q  to  Pi 
and  ddiQ)  is  the  angle  between  the  two  rays  bounding  a 
tube.  As  Pi  tends  to  Q,  s  tends  to  zero  and  Ai(Pi)  be- 
comes infinite  in  such  a  way  that  the  following  limit 
exists:  lim.,^o^tl,(/'i")  =  -'l/(0-  For  a  line  source  con- 
sisting of  a  unit  force  per  unit  length  in  the  .v  direction 
in  an  infinite  homogeneous  medium,  .1  /{Q)  is  given  by 
the  following  expression,'  in  which  we  write  /l/(.r)  to 
indicate  the  direction  of  the  force : 

^^x(0  =  e,./4(8^^)-5[r^(0]-J[p((3)]-isin&(()),     (15) 

.4/(0  =  e-/H8,ra;)-iCi'.(0]-Tp(0]-'cose((>)-    (16) 

In  Eqs.  (15)  and  (16),  e{Q)  is  the  angle  that  the  ray 
makes  with  the  positive  y  axis  at  Q.  The  unit  vector  n 
is  oriented  so  that  kX  n  points  in  the  positive  z  direction. 
In  accordance  with  the  geometrical  method,  we  assume 
that  Ai'^iQ)  is  also  given  by  Eqs.  (15)  or  (16)  for  in- 
homogeneous  and  bounded  media.  Now,  Eqs.  (4)  and 
(5)  yield  with  Pi  and  P  replaced  by  Q  and  Qa, 
respectively, 


r  dd{Q) 

4 


driQ.i)  piQ)p(Q,i)Vi{Qi,). 


sm 


e,m-    (17) 


Here,  and  in  subsequent  equations,  sin  is  to  be  used  for 
i=c  and  cos  for  i  =  s. 

Upon  inserting  Eq.  (17)  into  Eq.  (10),  we  obtain 


A,i^{P)==e''"K8Trio)-n 


-dr{Pi,)  p{P,j)  Vj{P.i) 


X- 


driP)     piP)     V,(P 

p'iQii)  VuiQij)  de{Q) 


p'{P.,)Vu(P„)dr(Q.i)  p{Q)p(Q 


- — T 


sm 


xlvmVRiiPij)E.(Qu)     6.AQ)-    (18) 


Here,  d^  is  the  angle  between  k,,-  and  the  y  axis,  and  we 
iiave  utilized  the  fact  that  Hi,  lies  in  the  xy  plane  and  is 
perpendicular  to  k,s.  Similarly,  we  find 

«.(0  =  -4cc.(O)C-sine,e(())>^"«3' 

+-4.c(Q)C-sin9,,(e)>'^-('3' 
+  .4..(e)[cose„(Q)>'>-(e' 

-f.l.,((2)[cos9..(0>--"''2).     (21) 

To  obtain  u,{Q),  we  use  Eqs.  (18)  and  (19)  with  P 
and  (2  interchanged,  taking  account  of  the  fact  that  the 
rays  from  Qio  P  are  also  the  rays  from  P  to  Q.  We  still 
let  i  denote  the  type  of  ray  from  Q,j  to  Q  and  j  denote 
the  type  of  ray  from  P  to  P,j  so  that  the  labeling  of  the 
points  Ptj  and  Qa  is  unchanged.  Then,  we  find  from 
Eq.  (19)  that  ^ij{P)=ifj,(Q).  Since  ii^Q)  results  from 
a  force  in  the  y  direction  at  P,  the  coefficient  Ai^iP) 
is  obtained  from  Eqs.  (15)  and  (16)  by  replacing  d  by 

We  now  insert  Eqs.  (20)  and  (21)  into  Eq.  (13). 
The  resulting  equation  must  hold  for  all  points  P  and  Q 
and  all  values  of  ic.  Therefore,  the  coefficients  of  e'«"j<C' 
ande''^"'^*  must  be  equal.  This  yields  the  four  equations 

AcciP)  cos9,,{P)=-A,c{Q)  sineUQ) ;        (22) 

A,,{P)  cos^.c {P)  =  A „ (Q)  cose„ (<3) ;  (23) 

AcsiP)  sine,s{P)=-A,c{Q)  sind^ciQ) ;        (24) 


A,,{P)smejP)  =  A.,{Q)cose,,{Q). 


(25) 


In  these  equations,  we  now  use  Eq.  (18)  and  the  cor- 
responding expression  with  P  and  Q  interchanged, 
noting  that,  in  interchanging  P  and  Q,  d,j{Q)  must  be  re- 
placed by  e,j{P)-{ir/2).  We  then  find  that  the  trigono- 
metric factors  and  the  minus  signs  in  these  equations 
cancel  out  and  that  they  all  become  of  the  form  Aij(P) 
=  Aji(Q),  with  the  sines  and  cosines  omitted.  We  now 
use  Eq.  (18)  for  A  ij{P)  and  for  Aji{Q),  equate  these  two 
expressions  omitting  the  trigonometric  factors,  and 
transpose  all  factors  involving  P  and  Pa  to  one  side 
and  those  involving  Q  and  Qa  to  the  other.  After  simpli- 
fying, we  can  write  the  resulting  equation  in  the  form 


Inserting  Eq.  (14)  into  Eq.  (11)  yiek 


/       ds/Vn+i.  /      ds/Vi.  (19)      L  I  </r(P)  )  Q I  de{P)  \  pA 

J  Qii  J  P,i  ,w„    ■ 


When  Eqs.  (18)  and  (19)  are  used  in  Eq.  (12),  they  yield 
u(P).  The  y  component  of  u(P)  is 

Hy{P)  =  AUP)cosecc{P)e"'"^''^ 

+^..(P)cos9.c(/')e,'^"'^' 
+/l„(P)sine.,(P)e"'-(0> 

+As,{P)smeUP)e'^"^''\     (20) 

'  H.  Lamb,  "On  the  Propagation  of  Tremors  over  the  Surface 
of  an  Elastic  SoUd,"  PhiL  Trans.  Royal  Soc.  (London)  203,  1-42 
(1904). 


driP^]     (rfr(P,y)] 
(P)  \q[  de{P)  ,P. 

VAP,MP^J)LViiP)yRiiP'i) 

Vr{P.,)p'{P,,)        E,iP„) 
driQ.i)]     \driQ,j)]    nJ 


=[ 


driQ)  J  pi  dd{Q)  IqJ 

■.iQ,MQu)LvmyR.iQ. 


X- 


VRiQu)p'(Q,i)       E,{Q,i) 


-.     (26) 


In  Eq.  (26),  the  subscript  Q  or  P  on  a  geometrical  ratio 
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indicates  the  source  of  the  rays  to  be  used  in  calcuhiting 
that  ratio. 

The  geometrical  ratios  in  Eq.  (26)  can  be  calculated 
simply  by  letting  Q  approach  Q,j  so  that  the  ray  from 
Q  to  Q,j  is  essentially  straight,  and  similarly  for  J'  and 
P.j.  Then,  {(/7(<2„)/(/e(())}Q  =  </((?,i,(?)  i^  t'lif  distance 
from  Q,j  to  Q,  which  will  be  complex  if  Q,;  is  complex. 
Furthermore,  from  Ref.  l,Eci.  (10),  we  find  when  the 
rays  are  straight  that 


1  driQ) 


=  [l  +  (K+e,)</(y,„0  sece,]^,,,-".     (27) 


Here,  k  denotes  the  curvature  of  llie  surface  at  1^',^,  9, 
is  the  value  of  0  given  by  Eq.  (1)  at  Q„  with  1'=  V.UJ.j), 
and  di  is  the  derivative  of  dt  with  respect  to  arclength 
along  the  surface.  I'Vom  Eq.  (27)  and  the  preceding 
ratio,  we  have 


driQ.j)]     \dr{Q,i 


=  [K-o+(K+e,)sece,]g.,-'.      (28) 


It{q)  \p[de(Q)  U 

Here,  ko=  l/d{Q„,{>}  is  the  curvature  of  the  wavefront 
incident  at  Qij. 

We  now-  use  Eq.  (28)  on  the  right  side  of  Eq.  (26), 
and  on  the  left  side  we  insert  the  corresponding  expres- 
sion with  P's  and  Q's  interchanged  and  with  dj  in  place 
of  di.  Then,  the  left  side  of  Eq.  (26)  is  found  to  be  in- 
dependent of  Q  and  the  right  side  independent  of  P. 
Therefore,  each  side  is  a  constant,  which  we  shall  write 
as  a~~.  Furthermore,  since  we  have  let  P  approach  P„, 
we  may  replace  P  by  P,j  on  the  left  side  and  similarly 
Q  by  Qij  on  the  right.  Upon  equating  the  left  side  of 
Eq.  (26)  to  a~',  we  obtain  an  equation  for  Ej{P,j)  that 
yields 

Ei{P,j)  =  a-[_Ko+  (K+Bi)  sec^Jp, -![!'; (/'„)]^ 

XRjiP,i)p(P.,)/VH(P,,)p'(I'u)-     (29) 

It  is  convenient  to  introduce  the  new  radiation  coef- 
ficient Rj'iP.j)  defined  by 

R/{P,j)  =  aR,iP.jVAPv)y 

XCp'(/^,)!•«(/^,)]-i[p(/'.;)]*•     (-^0) 

Then,  Eq.  (29)  becomes 
E,iP.j)  =  a[Ko+  (K+dj)  sece,]/>,,-5 


X- 


LVAPM 


Lv^iP.j)j 


■PJP.,) 
■p'iP.d-i 


(M) 


This  is  the  relationship  between  Ej  and  Rj  that  we  have 
been  seeking.  We  see  from  Eq.  (26)  that  it  also  holds 
when  Pij  is  replaced  by  Q,j. 

Although  Eq.  (31)  has  been  derived  in  the  special 
case  of  cylindrical  symmetry  (i.e.,  for  the  two-dimen- 
sional case),  we  assume  that  it  holds  in  general,  provided 
that  Ko  and  k  are  the  curx'atures  of  the  incident  wave- 


front  and  the  surface  in  the  plane  of  incidence.  This  is 
the  plane  containing  the  tangent  to  the  incident  ray  and 
the  normal  to  the  surface.  In  addition,  6j  is  the  deriva- 
tive of  dj  along  the  surface  ray. 

Let  us  now  use  Eq.  (31)  for  EliiQi,)  and  eliminate  £, 
from  Eq.  (10).  In  this  way,  we  obtain 

A„{P)  =  A.{Q„)Rj'{P„)Ri{Q„) 

■  r V.JQ.MQ.j)  dajQ.,)  dT{P„)-}i 

V,{P)p(P)    da{P„)  dT{P)\ 

■[Ko+(^  +  e,)sec6>,]o„-i.     (32) 


X 


rViiQ^ 


When  lv|.  {?<2)  is  used  in  Ec).  (12)  ,it  yields  the  surface- 
wave  field  at  P  due  to  an  arbitrary  wave  incident  upon 
the  surface.  We  note  that  Eq.  (32)  does  not  involve  p', 
which  has  been  absorbed  into  the  radiation  coefficient 
Ri .  In  the  next  section,  we  determine  /?,'  for  a  free 
surface,  using  a  method  that  can  also  be  applied  to  other 
kinds  of  surfaces.  Once  Ri  is  know-n,  Eqs.  (12)  and  (32) 
will  express  the  surface-wave  field  in  terms  of  known 
geometrical  quantities  and  material  properties. 

For  a  homogeneous  medium  bounded  by  a  plane 
surface,  the  space  rays  and  surface  rays  are  straight 
lines.  Then,  Eq.  (3)  yields,  with  v'.(<3)  =  0, 

V''.;(P)  =  C0rf((),C,;)/  \\  +  0.hHQ,j,P„)/  Vr 

+o.hI{P,„P)/\-j.     {.U) 

In  the  two-dimensional  case  of  a  homogeneous  medium, 
Eq.  (32)  becomes 

.4„(P)  =  /?//?/(F,/F,)S[.„(()./)]-M.(0,y).     (34) 

In  Eq.  (34),  we  have  made  use  of  the  facts  that  ^,=  0, 
K  =  0,  (h(Q,j)/daiP,,)=l,  dT{P,j)/dT{P)=\  from  Eq. 
(27),  and  all  coefficients  are  constants.  If  the  field  is 
produced  by  a  line  source  at  Q,  it  follows  from  Eqs.  (4) 
and  (5)  that  [ko((?.;)]"*-1,(<3.j)  does  not  depend  upon 
the  distance  from  Q  to  Qij  but  depends  only  upon  the 
direction  d,j{Q)  of  the  ray  from  Q  to  Q,j.  Calling  this 
product  .1  /[^,j((})],  we  may  write  Eq.  (34)  as 


[.,(/')  =  i?/A'/(F,/F,)U/[^„(())]. 


(35) 


The  quantity  A  ,'\J),j{Q)~\  expresses  the  strength  and 
angular  pattern  or  radiation  pattern  of  the  source  and 
depends  upon  the  source  characteristics.  In  the  case 
of  a  unit  force  per  unit  length  in  the  .t  direction,  A  '  is 
given  by  Eqs.  (15)  and  (16).  Upon  using  Ecis.  (21) 
and  (22)  in  Eq.  (12),  we  obtain  u(/')  due  to  a  line 
source  at  Q  in  a  semi-infinite  homogeneous  medium. 

IV.  DETERMINATION  OF  THE  RADUTION 
COEFFICIENT 

We  now  use  the  preceding  results  to  determine  the 
field  u(/')  due  to  a  compressional  line  source  at  Q  in  a 
semi-infinite  homogeneous  medium  bounded  by  a  free 
surface.  Then,  we  compare  the  results  with  the  exact 
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Fig.  2.  An  elastic  medium  with  a  free  surface  at  y  =  0  occupies 
the  half-space  >'>0.  The  compressional  Une  source  is  parallel  to 
the  z  a.xis  and  intersects  the  plane  :  =  0  at  Q  with  coordinates 
x  =  0,  3'  =  )'o.  There  are  two  ray  paths  from  the  source  to  the  ob- 
servation point  at  P  and  are  denoted  by  QQcPcP  and  QQcP.P. 


written  |.t|  instead  of  x,  so  the  results  will  be  valid  for 
x<Q.  Since  we  consider  a  compressional  source, 
,4  J  (Q)  =  0,  so  we  needn't  consider  the  shear  rays  from  Q. 
Let  us  now  use  Eqs.  (38)  and  (39),  together  with 
Eqs.  (35)  and  (12),  to  obtain  for  the  surface  wave  at  P 

u(P)  =  ^/(e.)i?c'=kcexp{wrrl|.rIsinec 

+  (yo+y)  cos^e]}+.-i/(ee)-Kc'i?/(rc/r.)4n, 

Xexp{;<jrr'CI  ^1  sin^c+yocos^c 

^-Vr'V.ycosd,-]}.     (40) 


The  unit  vectors  k,,  and  n,  are 

kc=  (sin^c,  cosdc), 
n,=  (cos^s,  —  sin^s). 


solution  of  the  equations  of  elasticity  for  this  case,  when 
that  solution  is  expanded  asymptotically  for  large  values 
of  the  frequency  co.  This  comparison  will  show  that  the 
two  results  are  identical  when  the  radiation  coefficients 
Ri  are  chosen  appropriately.  We  then  assume  that  these 
coefficients  have  the  values  so  determined  in  all  cases 
involving  a  free  surface,  provided  the  material  constants 
in  the  expressions  for  7?/  are  given  the  values  appropri- 
ate to  the  point  of  excitation  or  radiation. 

We  begin  by  observing  that  the  point  Q,j  depends 
only  upon  i  and  Q,  while  the  point  P,y  depends  only 
upon  J  and  P.  This  is  because  the  ray  from  Q  to  Q,j  is 
a  straight  line  that  makes  an  angle  with  the  normal  to 
the  surface  determined  by  Eq.  (1).  In  Eq.  (1),  only  Vi 
occurs,  so  the  angle  depends  on  F,,  and,  therefore,  the 
location  of  Qa  depends  only  upon  F,  and  the  location 
of  Q.  The  same  considerations  apply  to  P.j.  Therefore, 
we  write  Q,  instead  of  Q^j  and  Pj  instead  of  P.y. 

Let  the  medium  occupy  the  half-space  y>0  with  a  free 
surface  at  y=0.  Let  the  line  source  be  parallel  to  the  z 
axis  and  let  it  intersect  the  plane  z=  0  at  ()  with  coordi- 
nates x=Q,  y  =  yo.  Then,  the  coordinates  of  Qc  are 
(yo  tan^c,  0).  If  the  observation  point  P  is  at  {x,y)  with 
.T>0,  then  Pc  is  at  {x-yia.ndc,  0)  and  P.  is  at 
(.X-— y  tan9„  0)  (see  Fig.  2).  Between  these  points,  we 
find  the  distances  to  be 

d{Q,Qc)  =  yai>ecec\ 

d{Qc,Pc)  =  .V-  (vu+y)  tane, ;  (36) 

d{Pc,P)  =  y&ccec; 

diQc,P,)  =  x-y„  iandc-y  tan^, ;  (37) 

d{Ps,P)  =  yse^e,. 
Upon  using  Eqs.  (36)  and  (37)  in  {33),  we  obtain 
'Pec  (P)  =  CO  I  •,-'[  I  .V I  miec+  (yo+y)  co^Oc'] ;  (38) 

'Pcs{P)  =  oo\\-'l\x\sine,+  yocose,+  Vc-H\ycoses'].  (39) 

Here,  we  have  used  the  relations   V,=  VRsmd,  and 
yc=  I'«  sin^c,  which  follow  from  Eq.  (1).  We  have  also 


(41) 
(42) 


The  result  (40)  involves,  in  addition  to  known  geo- 
metrical quantities,  the  three  unknown  quantities 
A/{Oc),  Re',  and  R,'.  Let  us  consider  the  exact  solution 
for  the  displacement  in  an  infinite  homogeneous  medium 
due  to  a  compressional  line  source  of  unit  strength  (unit 
volume  per  unit  length).  This  solution,  when  expanded 
for  large  values  of  a',  yields 

A/{de)={W»^Vey-e"-"-  (43) 

Next,  we  consider  the  exact  solution  for  the  displacement 
due  to  the  same  source  in  a  semi-infinite  homogeneous 
medium  with  a  free  surface.  This  solution  has  been  ex- 
panded asymptotically  for  large  values  of  o).  In  the 
notation  of  the  previous  paragraphs  of  this  section,  the 
surface  wave  at  P  is  found  to  be 

viix,y)^Kio,/Ve)e"'H2sm^0e-Ve'V,--r-LdG(e,)/dey' 
Xexp{Ja)FrTi-'*-"|sine,+  (yo+y)cos6iJ} 
+n,(a)/F,)e"/22  sinec(2  sin-6i,-  Ve-Vr") 

X  {Ve-v,---^m-e,y>idG{d,)/de']-' 

Xexp{toFc"^[|.T|sinec+yocosec 

+  Fr'F.y  cose  J}.     (44) 

In  Eq.  (44),  the  function  G(d)  is  defined  by 

6(6)=  (2  ain-e-  F/Fr=)=+4  mrB 

Xcose(F/Fr--sin20)i.     (45) 

Upon  comparing  Eq.  (44)  with  Eq.  (40),  we  find  that 
they  agree  exactly,  provided  that 

A/{dc)Re"  =  e'"HWye){2sm'dc-VJ'Vr'y 

xldG{d,}/dd-y\    (46) 

A/{ee)R/R/iVc/V.)^=2e-'-{o^/y>) 

Xsin6i<:(2  sin-^c-  rc-F,--)(Fc=Fr--sin=6l,)' 

XldG(d,)/ddyK     (47) 

When  Eq.  (43)  is  used  for  .4  /  (e,)  in  Eqs.  (46)  and  (47), 
they  can  be  solved  for  P/  and  R/.  The  results  are 

R,'=  (87ra)/F,)ie— '8(2  sin^^.-  r,=  rr=) 

xLdGie.ydd'j-K   (48) 
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/?/=  (87ru)/Kc)'e~''"'*2  sin^c  cosOc  face  wave  excited  on  a  curved  free  surface  by  an  arbi- 

X  [,dG {dc)/d9']^^  ( Vc/  F.) '.     (49)      trary  incident  field  in  an  inhomogeneous  elastic  medium. 

It  is 
The  above  equations  determine  the  radiation  coefficients. 

Let  us  now  insert  Eqs.  (48),  (49),  and  (32)  into  Eq.      u(/^)  =  \/7re-'<'/^'a)i[B„(P)e"^^"''^'+B,c(P)«-'"^"'<^> 
(12)  andreplace  the  amplitudes  AdQcj)  a.nd  A,(Q,j)  by  +B,,(P)e'"«'"'(^)+B,.f'''"'(^'],     (50) 

Eqs.  (4)  and  (5).  This  yields  the  final  result  for  the  sur-      where 


r  I  '^  (Qcc)p  {Qcc)da  (Qcc)dr  (Fee)  -\ » 
B..(/')  =  2(2)M.(e..)[F.(/'..)]-i  /  •[2sin->e.-lVr-]p„ 

X[2  sin^Sc-  V^Vi^Qj .[k„+  (K+^c)]o„-*kc.(/') ;  (51) 

L    dd    Jp„L    dd    Jo„ 

r^%(e..)l"Yl'.((?»c)p((?,c)</a(e„e)<fr(P..)-lS 

B..(i')  =  2(2)tl.(C.c)CF.(Fj]H    — ^-  ^ -[2  sin=0.-r/rr=]p„ 

LK,(e.,c)J      L      Fc(/')p(P)'/c^(P.r)</r(P)     J 

,     raG(e,)-]-»  r6G(e.)-|-i 
X  [2  sine,  cose,]Q„    •  [k..+  {K+6.)-]Q.r^KAP) ;  (52) 

L   de  Ap,X  dd  Jq,c 

rVc{Pcs)-\"'rVciQc,)piQcs)da{QcMTiPc.)-]^ 

B.„(p)  =  2(2)M,(ft.)[i^(p„)]-i  — ^^ — ^ -[2  sin=e.- r.^rr^]^,. 

Lf,(P,,)J      L      VAP)p{P)dc{PcMr{P)     J 


X[2sine,cos6».]r,.    •[^o+('^+(',)]o„-*n„(P) ; 

L   dd   AqX   dd   J/.„ 

T>  ,r.s    ...^;.  .^  Nr,.  ..  s.  jyc{Pss)^'i'r\\{Qs,)-t"  r\\{Q.s)p{QsM<^{Q.,)dT{Ps,)-\^ 
B„(P)  =  2(2)M,((),.)[F.(P..)]-M — ■    

Lf,(p„)J    IvXQss)^     L    i',,(p)p(p)</<7(/',,)rfr(P)   J 


(53) 


■  [2  smffc  cos9,]p„[2  sin^,  cos^c].?,.-    

L   de  Jo.x   dd  . 


■['co+(K+0,O]o..-*n„„(/');  (54) 


and  the  phase  ipi/iP),  which  is  independent  of  6.',  is  It  is  possible  to  use  Eq.  (55)  in  its  present  form  to 

given  by  ipij{P)/ui.  find  the  waveform  for  transient  waves.  However,  it  is 

simpler  to  use  another  procedure  that  makes  use  of  the 

V.  TRANSIENT  WAVES  special  form  of  the  frequency  dependence  obtained  in 

The  results  obtained  in  the  preceding  section  apply  E_q.   (50).  The  time-harmonic  solution  given  by  Eq. 

to  time-harmonic  waves.  We  now  transform  them  to  ^     '  '^       ^      form 

apply  to  transient  waves.  Suppose  that  we  have  a  time-  u(x,/)  =  \/7re~'*'''^'co'  ^  B.yc'"*'''',            (57) 

harmonic  solution  u(x,a))e~'"'  of  the  elastic  equations  *.J 

of  motion  for  a  source  distribution  with  harmonic  time  ^here  the  summation  consists  of  four  terms  j=c,  5; 

dependence.  Then  it  is  easily  shown,  by  means  of  Fourier  j^^^  s.  The  amplitudes  B„-  and  phases  ^,/,  which  are 

transforms,  that  the  time-dependent  solution  U(x,/)  independent  of  frequency,  are  easily  found  by  compari- 

for  the  same  source  distribution  with  an  arbitrary  time  ^^n  of  Eq.  (57)  with  (50)-(54).  Using  Eq.  (57)  in  (55), 

dependence  F{t)  is  given  by  ^e  obtain 


U(x 


,0^/'"u(x,co)/(.).-...  (55)      U(x,0  =  x/.e-/^.i:B,/_^.^.—.n/(.)...     (58) 

,,         ,,,.,„.  ,  TT^/N       ,  Now,  we  differentiate  Eq.  (56)  a  times  with  respect  to  / 

Here,  ]{w)  is  the  Fourier  transform  of  F{t)  and  H    bt  ' 

J{i^)e'''"'du.  (56)  —FiO^i-i)"        co''e-'"'/(co)(fco.  (59) 

dt"  J-^ 
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If  we  letQ:  =  5,  then  a  comparison  of  Eq.  (58)  with  (59) 

shows  that  ^"'^                                ^       ,, 

U(x,/)  =  n/tt  Z  liii{d/dt)iF{l-  .pi/).          (60)  U(x,/)  =  Z  Bi;—  /    it*-T)-iF{T)dT,        (64) 

i,i  i.j       at*  J  0 

Thus,  the  desired  transient  waveform  for  a  given  input  ^^j^^^.^  ^*  ^  ^_  _^ . .,  j^  ^^^  delayed  time.  The  above  equa- 

Fit)    is    given    by    the    one-half-order    derivative    of  tions  yield  the  transient  response  for  an  arbitrary  pulse 

F{t->pi/).  It   is  well-known   that   the  one-half-order  /?(;)  ^nd  are  the  main  results  of  this  section, 

derivative    is    given    by    either    of     the    following  ^^  ^  ^^^^^^  example,  suppose  that  F{l)  is  the  unit 

expressions^:  step  function : 

(0,    /<0 

, . .                       .                        (65) 

V^rJo                dr  ll,     />0 


(-)F{t)  =  —  I   {l-T)-^-F{r)dT.  (61)  F{t)=\ 

\dU  \/iTjo  d-r  i- 

/(Z\i  1     d    r< 

(-)F{t)  = /    {t-r)-^F{T)dr. 

\di/  x/wdtJa 


•d\^             Id    f'  Then,  (/F(/)/a'/  =  6(/),  and  Eq.  (63)  yields  immediately 

"'"   = /    {i-T)-^F{T)dT.         (62) 

\/^  dt  Jo  U(x,0  =  Z  Bi;0-  ^,/)-i.                  (66) 

Upon  substituting  Eq.  (61)  or  (62)  into  (60),  we  obtain 

two  different  expressions  for  the  transient  waveform.  Other  pulses  may  be  treatefl  in  a  similar  manner. 
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Abstract 

The  amplitudes  of  the  fields  in  the  shadows  of  two  screens  with  rounded 
ends  are  deternaned  as  functions  of  bA.  Hei^  b  denotes  the  radivis  of  curvature 
of  the  end  of  the  screen  and  X  denotes  the  wavelength  of  the  normally  incident 
plane  wave.  For  one  of  the  screens  the  amplitude  decreases  monotonically  as 
bA  increases  irtien  the  field  vanishes  on  the  screen.  But  if  the  normal 
derivative  vanishes  the  amplitude  first  oscillates  as  bA  increases  and  then 
finally  decreases.  For  the  other  screen  -  the  more  realistic  one  -  the 
amplitude  probably  decreases  monotonically  in  both  cases  as  bA  increases. 

The  calculations  are  based  on  the  geometrical  -Uieory  of  diffraction.  For 
the  second  screen  a  calculation  from  the  exact  series  solution  is  also  made. 
The  close  agreement  between  the  two  results  is  a  confirmation  of  the  geometrical 
diffraction  theory. 
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1,  Introduction 

Because  of  diffraction,  there  is  always  some  illumination  in  a  shadow. 
The  intensity  of  this  illumination  depends  upon  the  wavelength  X  of  the 
incident  radiation  and  upon  the  size  and  shape  of  the  shadow-forming  object. 
In  order  to  examine  this  dependence,  we  have  considered  the  shadows  cast  by 
two  semi-infinite  screens  with  rounded  tips  of  radius  b  upon  either  of  which  a 
plane  wave  of  unit  amplitude  is  normally  incident.  (See  Figures  1  airi  2).  We 
find  that  in  either  case,  at  a  large  distance  r  from  the  tip  the  diffracted 
field  u  in  the  shadow  has  the  form 

Kkr*") 
(1)  u  -  5 IL.  f(e,kb)  . 

In  (1)  k  ■  2n/\  is  the  propagation  constant  and  f(0,kb)  is  the  diffraction 
pattern  or  complex  amplitude  of  the  cylindrical  wave  scattered  from  the  tip  in 
the  direction  0.  The  angle  0  is  measured  from  the  shadow  boundary. 

In  Figure  3j  the  amplitude  |f(0,kb)|  is  shown  as  a  function  of  kb  for 
0  -  n/k  for  the  screen  in  Figure  2,  The  encircled  points  show  that  the  shadow 
cast  by  this  screen  becomffl  darker  as  kb  increases  when  u  ■  0  on  the  screen. 

Thus,  the  shadow  becomes  darker  as  the  radius  of  the  tip  increases  as  compared  to 

-1/2 
the  wavelength.  The  factor  k  '  in  the  denominator  of  (1)  alao  makes  the 

shadow  become  darker  as  the  wavelength  decreases,  independently  of  the  radius  of 

the  tip.  On  the  other  hand,  if  du/3n  »  0  then  |f(T-,kb)|  oscillates  as  kb 

increases.  Ultimately  it  ceases  oscillating  and  decays  monotonically  as  (6) 

shows.  This  oscillation  is  due  to  interference  between  the  two  diffracted  rays 

which  occur  in  this  case.  Although  the  same  two  rays  occur  when  u  ■  0  on  the 

sci-een,  one  of  them  decays  so  rapidly  that  it  produces  negligible  interference 

and  no  oscillation. 
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Figure  1  Cross-section  of  a  screen  of  width  2b  with  a  rounded  tip  of  radius  b. 
A  plane  wave  is  normally  incident  upon  it  from  the  left.  The  dotted  line  is 
the  shadow  boundaiy.  An  incident  ray  which  grazes  the  tip  is  shown  together 
with  one  of  the  diffracted  rays  iriiich  it  produces  in  the  shadow.  The  angle  between 
the  diffracted  ray  and  the  shadow  boundary  is  0. 


Figure  2  Cross-section  of  an  infinitely  thin  screen  with  a  cylindrical  tip  of 
radius  b,  A  plane  wave  is  normally  incident  upon  it  from  the  left.  The  dotted 
line  is  the  shadow  boundary.  An  incident  ray  which  grazes  the  tip  is  shown 
together  with  the  two  diffracted  rays  which  it  produces  making  the  angle  6 
with  the  shadow  boundary. 
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An  important  feature  of  Figure  3  is  the  close  agreement  between  the 
predictions  of  the  geometrical  theory  of  diffraction  and  the  exact  solutions  of 
the  boundary  value  problems  for  the  field. 

In  Figure  h  the  amplitude  |f(©,kb)|  is  shown  as  a  function  of  kb  for 
e  »  n/k  for  the  screen  of  Figure  1.  From  the  physical  significance  of  the  approxi- 
mations on  which  this  figure  is  based,  it  is  likely  that  only  the  decreasing 
parts  of  the  curves  are  correct.  For  this  screen  the  amplitude  probably  decreases 
monotonically  for  all  values  of  kb. 

2.  Analysis 

If  the  radius  of  the  tip  is  zero  either  screen  is  a  half -plana.     In 
this  case  Sommerfeld's  well-kiK)wn  expression  for  u  is  of  the  form  (1)  with  f 
given  by 

(2)  f(©,0) 


sec  "s(9+n)  t  esc 


^(e*n)j  . 


2  v/H 

The  upper  sign  in  (2)  applies  to  a  scr«en  on  which  u  «  0,  while  the  lower  sign 
applies  to  one  on  which  8u/8n  ■  0.  At  the  other  extreme,  when  kb  is  large,  the 
field  in  the  shadow  can  be  determined  by  means  of  the  geometrical  theory  of 
diffraction''  -^ .  In  Figures  1  and  2  are  shown  the  one  or  two  diffracted  rays 
which  go  in  a  given  direction.  By  the  methods  and  formulas  in  [ij  we  find 
that  far  from  the  tip  the  field  associated  with  the  ray  in  Figure  1  is  given  by 


^  .  'J^    .(..,1/3  !»  S'""'"'"'  e^(-  ♦  fe  -  r' 


(3) 


This  same  expression  applies  to  the  upper  ray  in  Figure  2,   while  for  the  lower 
ray  we  find 

(U)         u,  -  »  5 1.  .  (kb)^/5  -2  e  °  ^i(kb|n-e|  ♦  T5  -  ^  )  . 
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The  upper  sign  in  (U)  applies  if  u  -  0  on  the  screen,  while  the  lower  sign 
applies  if  8u/6n  ■  0,     In  these  two  cases  the  constants  C     andT    are  given  by 

u-0  |H.o 

C  .910719  1.53228 

o 

-r^      1.855757  e^"/^  .308617  e^"/^    . 

The  fields  u  ■  tl  and  u  -  n,  ♦  u^  are  of  the  form  (l)  with  f  given  by  the 
respective  formulas: 

^irkb+r^(kb)^/3l€ 

v/2 


(5)  f(e,kb)  .  (kb)^/3   ^  e-i"/6  ,i|kb.rjkb)-^^|e 


(6)  f(9,kb)  -  (kb)^/3    -2  e-^"/*^ 

1/2 


i 

e 


[kb+T^Ckb)^/^]  ;^i|kb+r^(kb)^/^](n-e) 


In  Figures  3  and  h  the  values  of  |f(T-,kb)|  are  shown  based  on  (2)  for  kb  -  0 
and  on  (5)  and  (6)  (solid  line)  for  kb  >  .1  for  screens  on  which  either 
u  «  0  (lower  graph)  or  3u/Sn  «  0  (upper  graph). 

For  the  screen  of  Figure  2  the  field  u  has  been  determined  exactly  for  all 

r2i 

values  of  kb  by  the  method  of  separation  of  variables  ^  -' .     The   result  for  f , 

which  coincides  with  (2)  for  kb  ■  0,  has  been  expanded  asymptotically  for 

kb  large  [3].  It  then  agrees  with  (6),  We  have  used  this  solution  to 

calculate  |fl  in  the  range  0  <  kb  <  10,  Ihe  results  are  shown  as  circled  points 

in  Figure  3  for  both  the  case  in  which  u  «  0  (lower  graph)  and  that  in  which 

du/3n  =  0  (upper  graph)  on  the  screen.  The  broken  line  in  the  upper  graph 

of  Figure  3  for  the  case  3u/3n  ■  0  is  based  upon  (6)  with  the  factor  kb  +?l(kb) 

in  the  exponents  replaced  by  the  more  accurate  expression  kb  ■♦■?T(kb)  '  -<r  (kb) 

The  constant  (T  has  the  value  L  J  <r  •  ,1U5U63  e"  "'•^. 
o  o 

The  agreement  between  the  circled  points  based  upon  the  series  solution 

and  the  simple  results  (6)  of  the  geometrical  theory  of  diffraction  ^  -^  is 
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another  confirmation  of  the  geometrical  theory.  The  labor  involved  in  computing 
the  curves  from  (6)  was  almost  negligible.  However,  the  calculations  of  the 
points  from  the  series  solution  required  a  great  deal  of  time.  Four  series  of 
Bessel  functions  had  to  be  computed  for  each  point  and  each  required  about 
1»2  kb  ♦  2  terms  to  jdeld  three  decimal  accuracy.  To  use  the  series  for 
kb  >  10  would  require  tables  of  Bessel  functions  of  high  order  as  well  as 
many  terms.  In  view  of  the  agreement  of  the  series  result  with  (6),  this  does 
not  seem  to  be  necessary.  Thus,  for  example,  we  may  conclude  from  (6)  that 
for  kb  large  |f |  ultimately  decreases. 
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298  Shames  Drive 
Brush  Hollow  Industrial  Park 
Westbury,  New  York 

Electromagnetic  Research  Corroratlon 
711  Ibth  Street,  N.  W. 

Washington  5,  D.  C, 
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Lockheed  Aircraft  Corporation 

2?5S  N.  Hollywood  Way 

California  Division  Engineering  Library 

Department  72-75,  Plant  A-1,  Bldg.  63-1 

Burbanlt,  California 

Attn:   N.  C.  Harnois 

The  Martin  Company 

P.  0.  Box  179 

Denver  1,   Colorado 

Attn  I     Mr.   Jack  McCormick 

The  Glenn  L.  Martin  Company 
Baltimore  3,   MarylatiJ 
Attn  I     Engineering  Library 
Antenna  Design  Group 

Maryland  Electronic  Manufacturing  Corp. 
?009  Calvert  Road 
College  Park,   Maryland 
Attn!     Mr.   H.  Warren  Cooper 

Mathematical  Reviews 
190  Hope  Street 
Providence  6,   Rhode  Island 

The  W.  L.  Maxson  Corporation 

li60  West  3lith  Street 

New  York,   N.  T. 

Jttn:     Miss  Dorothy  Clark 

McDonnell  Aircraft  Corporation 

Lambert  Saint-Louis  Municipal  Airport 

Box  516,   St.  Louis  3,   Missouri 

Attn:     R.  D.  DAtrlch,   Engineering  Library 

McMillan  Laboratory,    Inc. 

Brownville  Avenue 

Ipswich,   Massachusetts 

Attn:     Security  Officer,   Document  Room 

Melpar,    Inc.  (2) 

3000  Arlington  Boulevard 

Falls  Church,    Virginia 

Attn:     Engineering  Technical  Library 

Microwave  Development  Laboratory 
90  Broad  Street 
Babson  Park  57,   Massachusetts 
Attn:     N.  Tucker,   General  Manager 

Microwave  Radiation  Company  Inc. 

19223  South  Hamilton  Street 

Gardens,   California 

Attn:     Mr.  Morris  J.  Ehrllch,   President 

Chance  Vought  Aircraft,    Inc. 
931lt  West  Jefferson  Street 
Dallas,   Texas 
Attn:     Mr.  H.  S.  White,   Librarian 

Northrop  Aircraft,    Inc. 
Hawthorne,   California 

Attn:     Mr.  E.   A.  Freitao,   Library  Deptjll:5 
inOl  E.  Broadway 

Remington  Rand  Univ.  -  Division  of  Sperry 

Rand  Corporation 
1900  West  Allegheny  Avenue 
Philadelphia  29,    Pennsylvania 
Attn;     Mr.    John  F.   McCarthy 

R  and  D  Sales  and  Contracts 

North  American  Aviation,    Inc. 
1221U  Lakewood  Boulevard 
Downey,   California 
Attn:     Engineering  Library  U95-11S 

North  American  Aviation,   Inc. 
Los  Angeles  International  Airport 
Los  Angeles  1:5,   California 
Attn:     Engineering  Technical  File 

Page  Communications  Engineers,    Inc. 
710  Foxirteenth  Street,   Northwest 
Washington  5,   D.  C. 
Attn:     Librarian 

Philco  Corporation  Research  Division 

Branch  Library 

li700  Wlssachlckon  Avenue 

Philadelphia  hh,   Pa. 

Attn:     Mrs.  Dorothy  S.  Collins 


Plckard  and  Bums,    Inc. 

2b0  Highland  Avenue 

Needham  9U,   Mass. 

Attn:     Dr.   J.  T.  DeBettencourt 

Polytechnic  Research  and  Development 

Company,    Inc. 

202  Tlllary  Street 

Brooklyn  1,  New  York 

Attn:  Technical  Library 

Radiation  Engineering  Laboratory 
Main  Street 
Maynard,  Mass. 
Attn:  Dr.  John  Ruze 

Radiation,  Inc. 

P.  0.  Drawer  37 

Melbourne,  Florida 

Attn:  Technical  Library,  Mr.  M.L.  Cox 

Radio  Corp.  of  America 

RCA  Laboratories 

Rocky  Point,  New  York 

Attn:  P.  S.  Carter,  Lab.  Library 

RCA  Laboratories 
David  Sarnoff  Research  Center 
Princeton,  New  Jersey 
Attn:  Miss  Fern  Cloak,  Librarian 
Research  Library 

Radio  Corporation  of  America 
Defense  Electror^ic  Products 
Building  10,  Floor  7 
Camden  2,  New  Jersey 
Attn:  Mr.  Harold  J.  Schrader 

Staff  Engineer,  Organization 

of  Chief  Technical 

Administrator 

The  Ramo-Wooldrldge  Corporation 
P.O.  Box  U5U53  Airport  Station 
Los  Angeles  1:5,  California 
Attn;  Margaret  C.  Whltnah, 
Chief  Librarian 

Hoover  Microwave  Co. 
9592  Baltimore  Avenue 
College  Park,  Marvland 

Director,  DSAF  Project  RAND 

Via:  Air  Force  Liaison  Office 

The  Rand  Corporation 

1700  Main  Street 

Santa  Monica,  California 

Rantec  Corporation 

Calabasas,  California 

Attn:  Grace  Keener,  Office  Manager 

Raytheon  Manufacturing  Company 
Missile  Systems  Division 
Bedford,  Mass. 
Attn:  Mr.  Irving  Goldstein 

Raytheon  Manufacturing  Company 
Wayland  Laboratory,  State  Road 
Wayland,  Mass. 
Attn:  Mr.  Robert  Borts 

Raytheon  Manufacturing  Company 
Wayland  Laboratory 
Wayland,  Mass, 

Attn:  Miss  Alice  G.  Anderson, 
Librarian 

Republic  Aviation  Corporation 
Farmingdale,  Long  Island,  N.  Y, 
Attn;  Engineering  Library 

Thru;  Air  Force  Plant  Representative 
Republic  Aviation  Corp. 
Farmingdale,  Long  Island,  N.Y. 

Rheem  Manufacturing  Company 
9236  East  Hall  Road 
Downey,  California 
Attn:  J.  C,  Joerger 

Trans-Tech,  Inc. 
P.  0.  Box  3li6 
Frederick,  Maryland 


Ryan  Aeronautical  Company 
Lindbergh  Field 
San  Diego  12,  California 
Attn;  Library  -  unclassified 

Sage  Laboratories 
159  Linden  Street 
Wellesley  81,  Mass. 

Sanders  Associates 
95  Canal  Street 
Nashua,  New  Hampshire 
Attn:  N,  R.  Wild,  Library 

Sandia  Corporation,  Sandla  Base 

P.O.  Box  5R00,  Albuquerque,  New  Mexico 

Attn:  Classified  Document  Division 

Sperry  Gyroacope  Company 

Great  Neck,  Long  Island,  New  York 

Attn:  Florence  W.  Tumbull,  Engr.  Librarian 

Stanford  Research  Institute 

Menlo  Park,  California 

Attn;     Library,   Engineering  Division 

Sylvania  Electric  Products,    Inc. 
100  First  Avenue 
Waltham  5(i,   Mass. 

Attn;     Charles  A.  Thornhill,   Report  Librarian 
Waltham  Laboratories  Library 

Systems  Laboratories  Corporation 
1U^52   Ventura  Boulevard 
Sherman  Oaks,   California 
Attn:      Donald  L.  Margerum 

TRG,    Inc. 

17  Union  Square  West 

New  York  3,   N.  Y. 

Attn;     M.  L.  Henderson,   Librarian 

A.  S.  Thomas,  Inc. 

161  Devonshire  Street 

Boston  10,  Mass. 

Attn;  A.  S.  Thomas,  President 

Bell  Telephone  Laboratories 
Hurray  Hill 
New  Jersey 

Chu  Associates 
P.  0.  Box  387 
Whitcorab  Avenue 
Littleton,  Mass. 


Microwave  Associates, 
Burlington,  Mass. 
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Raytheon  Manufacturing  Company 
Missile  Division 
Hartwell  Road 
Bedford,  Mass. 

Radio  Corporation  of  America 
Aviation  Systems  Laboratory 
225  Crescent  Street 
Waltham,  Mass. 

Lockheed  Aircraft  Corporation 
Missile  Systems  Division  Research  Library 
Box  50li,  Sunnyvale,  California 
Attn:  Miss  Eva  Lou  Robertson, 
Chief  Librarian 

The  Rand  Corporation 
1700  Main  Street 
Santa  Monica,  California 
Attn:  Dr.  W.  C.  Hoffman 

Commander 

AF  Office  of  Scientific  Research 

Air  Research  and  Development  Command 

lUth  Street  and  Constitution  Avenue 

Washington,   D.   C, 

Attn:     Mr.  Ottlng,   SRY 

Westlnghouse  Electric  Corp. 
Electronics  Division 
Friendship  Int'l  Airport  Box  7I46 
Baltimore  3,   Maryland 
Attn;     Engineering  Library 


Wheeler  Laboratories,  Inc. 
122  Cutter  Mill  Road 
Oreat  Neck,  New  York 
Attn:  Mr.  Harold  A.  Wheeler 

Zenith  Plastics  Co. 
Box  91 

Gardena,  California 
Attn:  Mr.  S.  S.  Oleesky 

Library  Geophysical  Institute 
of  the  Dniverslty  of  Alaska 
College 
Alaska 

University  of  California 

Berkeley  h,   California 

Attni  Dr.  Samuel  Silver, 

Prof.  Engineering  Science 
Division  of  Elec.  Eng.  Electronics 
Research  Lab. 

nnlversity  of  California 
Electronics  Research  Lab. 
332  Cory  Hall 
Berkeley  U,  California 
Attn:  J.  R.  Whinnery 

California  Institute  of  Technology 
Jet  Propulsion  Laboratory 
bflOO  Oak  Grove  Drive 
Pasadena,  California 
Attn:  Mr.  I.  E.  Hewlan 

California  Institute  of  Technology 
1201  E.  California  Street 
Pasadena,  California 
Attn:  Dr.  C.  Papas 

Carnegie  Institute  of  Technology, 
Schenley  Park 

Pittsburgh  13,  Pennsylvania 
Attn:  Prof.  A.  E.  Helns 

Cornell  University 

School  of  Electrical  Engineering 

Ithaca,  New  York 

Attn:  Prof.  G.  C.  Dalman 

University  of  Florida 

Department  of  Electrical  Engineering 

Gainesville,  Florida 

Attn:  Prof.  M.  H.  Latour,  Library 

Library 

Georgia  Institute  of  Technology 

Engineering  Experiment  Station 

Atlanta,  Georgia 

Attn:  Mrs.  J.H.  Crosland,  Librarian 

Harvard  University 

Technical  Reports  Collection 

Gordon  McKay  Library,  303A  Pierce  Hall 

Oxford  Street,  Cambridge  3",  Mass. 

Attn:  Mrs.  E.L.  Hufschmidt,  Librarian 

Harvard  College  Observatory 
60  Garden  Street 
Cambridge  39,  Mass. 
Attn:  Dr.  Fred  L.  Whipple 

University  of  Illinois 
Documents  Division  Library 
Urbana,  Illinois 

University  of  Illinois 
College  of  Engineering 
Urbana,  Illinois 

Attn:  Dr.  P.  E,  Moyes,  Department  of 
Electrical  Engineering 

The- Johns  Hopkins  University 
Homewood  Campus 
Department  of  Physics 
Baltimore  U,  Maryland 
Attn:  Dr.  Donald  E.  Kerr 

Sandia  Corporation 
Attn:  Organization  lli23 
Sandia  Base 
Albuquerque,  New  Mexico 


Applied  Physics  Laboratory 
The  Johfls Hopkins  University 
8621  Georgia  Avenue 
Silver  Spring,  Maryland 
Attn:  Mr.  George  L.  Seielstad 

Massachusetts  Institute  of  Technology 

Research  Laboratory  of  Electronics 

Room  203-221 

Cambridge  39,  Massachusetts 

Attn:   John  H.  Hewitt 

Massachusetts  Institute  of  Technology 

Lincoln  Laboratory 

P.  0.  Box  73 

Lexington  73,  Mass. 

Attn:  Doonnent  Room  A-229 

University  of  Michigan 
Electronic  Defense  Group 
Engineering  Research  Institute 
Ann  Arbor,  Michigan 
Attn:   J.  A,  Boyd,  Supervisor 

University  of  Michigan 
Engineering  Research  Institute 
Radiation  Laboratory 
Attn:  Prof.  K.  M.  Siegel 
912  N.  Main  St., 
Ann  Arbor,  Michigan 

University  of  Michigan 
Engineering  Research  Institute 
Willow  Run  Laboratories 
Willow  Run  Airport 
Ypsilantl,  Michigan 
Attn:  Librarian 

University  of  Minnesota 

Minneapolis  lb,  Minnesota 

Attn;     Mr.   Robert  H.   Stumra,    Library 

Northwestern  University 
Microwave  Laboratories 
Evanston,    Illinois 
Attn:     R.  E.  Bean 

Ohio  State  University  Research  Found. 
Ohio  State  University 
Columbus  10,   Ohio 
Attn:     Dr.  T.E.  Tlce 

Dept.  of  Elec.  Engineering 

The  UnlTarsity  of  Oklahoma 
Research  Institute 
Norman,   Oklahoma 

Attn:  Prof.  C.  L.  F^rrar,  Chairman 
Electrical  Engineering 

Polytechnic  Institute  of  Brooklyn 

Microwave  Research  Institute 

5?  Johnson  Street 

Brooklyn,  New  York 

Attn:  Dr.  Arthur  A.  Ollner 

Polytechnic  Institute  of  Brooklyn 
Microwave  Research  Institute 
??  Johnson  Street 
Brooklyn,  New  York 
Attn:  Mr.  A.  E,  Laemmel 

Syracuse  University  Research  Institute 
Collendale  Campus 
Syracuse  10,  New  York 
Attn:  Dr.  C.  S.  Grove,  Jr. 

Director  of  Engineering  Research 

The  University  of  Texas 
Elec.  Engineering  Research  Laboratory 
P.  0.  Box  8026,  University  Station 
Austin  12,  Texas 
Attn:  Mr.  John  R.  Gerhardt 
Assistant  Director 

The  University  of  Texas 

Defense  Research  Laboratory 

Austin,  Texas 

Attn:  Claude  W.  Horton,  Physics  Library 

University  of  Toronto 

Department  of  Electrical  Engineering 

Toronto,  Canada 

Attn:  Prof.  G.  Sinclair 


Lowell  Technological  Institute 
Research  Foundation 
P.  0.  Box  709,  Lowell,  Mass. 
Attn:  Dr.  Charles  R.  Mingins 

University  of  Washington 

Department  of  Electrical  Engineering 

Seattle  5,  Washington 

Attn:  G.  Held,  Associate  Professor 

Stanford  University 
Stanford,   California 
Attn:     Dr.   Chodorow 

Microwave  Laboratory 

Physical  Science  Laboratory 

New  Mexico  College  of  Agriculture 

and  Mechanic  Arts 

State  College,  New  Mexico 

Attn:  Mr.  H.  W.  Haas 

Brown  University 
Department  of  Electrical  Engineering 
Providence,  Rhode  Island 
Attn:  Dr.  C.  M.  Angulo 

Case  Institute  of  Technology 

Cleveland,  Ohio 

Attn:  Prof.  S.  Seeley 

Columbia  University 

Department  of  Electrical  Engineering 

Momingside  Heights 

New  York,  N.  Y. 

Attn:   Dr.  Schlesinger 

HcGlll  University 

Montreal,  Canada 

Attn:  Prof.  0.  A.  Woonton 

Director,  The  Eaton  Electronics 

Research  Lab. 

Purdue  University 

Department  of  Electrical  Engineering 

Lafayette,  Indiana 

Attn:  Dr.  Schulta 

The  Pennsylvania  State  University 
Department  of  Electrical  Engineering 
University  Park,  Pennsylvania 

University  of  Pennsylvania 

Institute  of  Cooperative  Research 

3I4OO  Walnut  Street 

Philadelphia,  Pennsylvania 

Attn:  Dept.  of  Electrical  Engineering 

University  of  Tennessee 
Ferris  Hall 
W,  Cumberland  Avenue 
Knoxvllle  16,  Tennessee 

University  of  Wisconsin 
Department  of  Electrical  Engineering 
Madison,  Wisconsin 
Attn:  Dr.  Scheibe 

University  of  Seattle 

Department  of  Electrical  Engineering 

Seattle,  Washington 

Attn:  Dr.  D.  K.  Reynolds 

Wayne  University 

Detroit,  Michigan 

Attn:  Prof.  A.  F.  Stevenson 

Electronics  Research  Laboratory 
Illinois  Institute  of  Technology 
3300  So.  Federal  Street 
Chicago  16,  Illinois 
Attn:  Dr.  Lester  C.  Peach 
Research  Engineer 

Advisory  Group  on  Electronic  Parts 

Room  103 

Moore  School  Building 

200  South  33rd  Street 

Philadelphia  li,  Pennsylvania 
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Ionosphere  Research  Laboratory 

Pennsylvania  State  College 

State  College,  Pennsylvania 

^TTNi   Professor  A,  H.  Waynlcic,  Director 


Institute  of  Mathematical  Sciences 
2*^  Waverly  Place 
New  York  3,  New  York 
ATTN!   Librarian 

Electronics  Division 
Rand  Corporation 
1700  Main  Street 
Santa  Monica,  California 
ATTN:  Dr.  Robert  Kalaba 

National  Bureau  of  Standards 

Washington,  D.  C, 

ATTNi  Dr.  W.  K.  Saunders 

Applied  Mathematics  and  Statistir's  Lab. 

Stanford  University 

Stanford,  ailfomia 

ATTN:  Dr.  Albert  H.  Bowker 

Departmerit  of  Physics  and  Astronomy 
Michigan  State  College 
East  Lansing,  Michigan 
ATTN!   Dr.  A.  Leitner 

University  of  Tennessee 
Knoxvllle,  Tennessee 
ATTN:  Dr.  Fred  A.  Ficken 

Lebanon  Valley  College 
Annvllle,  Pennsylvania 
ATTN!  Professor  B.H.  Blsslnger 

General  Atomic 
P.  0.  Box  608 
San  Die, 
Attn:   1 


(3) 


12,   California 
Edward  Oerjuoy 


Department  of  Physics 
Amherst  College 
Amherst,   Mass, 
ATTNi     Dr.   Arnold  Arons 

California  Institute  of     Technology 
1201  E.  California  Street 
Pasadena,   California 
ATTN:     Dr.  A.  Erdelyi 

Mathematics  Department 

Stanford  University 
Stanford,  California 
ATTN!     Dr.    Harold  Levine 

University  of  Minnesota 

Minnpapolis  I'l,    Minnesota 

ATTN!      Professor  Paul  C.   Rosenbloom 

Department  of  Mathematics 

Stanford  University 

Stanford,   California 

ATTN:      Professor  Bernard  Epstein 

Applied  Physics  Laboratory 
The  Johns  Hopkins  University 
8621  Georgia  Avenue 
Silver  Spring.   Maryland 
ATTN:     Dr.  B.  S.  Gourary 

(2)Exchange  and  Gift  Division 
The  Library  of  Congress 
Washington  25,   D.  C. 

Electrical  Engineering  Department 
Massachusetts  Institute  of  Technology 
Cambridge  39,   Mass. 
ATTN:     Dr.  L.   J.  Chu 

Nuclear  Development  Associates,    Inc. 

5  New  Street 

White  Plains,  New  York 

ATTN:  Library 

California  Institute  of  Technology 
Electrical  Engineering 
Pasadena,  California 
ATTN:  Dr.  Zohrab  A.  Kaprielian 


Dr.  Rodman  Doll 
311  W.  Cross  Street 
Ypsllanti,  Michigan 

California  Inst,  of  Technology 
Pasadena,  California 
ATTN:  Mr.  Calvin  Wilcox 

Mr.  Robert  Brockhurst 

Woods  Hole  Oceanographic  Institute 

Woods  Hole,  Mass. 

National  Bureau  of  Standards 
Boulder,  Colorado 
ATTN:   Dr.  R.  Gallet 

Dr.  Solomon  L.  Schwebel 

3689  Louis  Road 

Palo  Alto,  California 

University  of  Minnesota 
The  University  of  Library 
Minneapolis  Ih,  Minnesota 
ATTN!  Exchange  Division 

Department  of  Mathematics 
University  of  California 
Berkeley,  California 
ATTN:   Professor  Bernard  Friedman 

Lincoln  Laboratory 

Massachusetts  Institute  of  Technology 

P.  0.  Box  73 

Lexington  73,   Massachusetts 

ATTN!     Dr.   Shou  Chin  Wang,   Room  C-351 

Melpar,  Inc., 

3000  Arlington  Boulevard 

Falls  Church,  Virginia 

ATTN:  Mr.  K.  S.  Kelleher,  Section  Head 

Hq.  Air  Force  Cambridge  Research  Center 

Laurence  G.  Hanscom  Field 

Bedford,  Mass. 

ATTN:  Mr.  Francis  J.  Zucker,  CRHD 

Hq.  Air  Force  Cambridge  Research  Center 

Laurence  G.  Hanscom  Field 

Bedford,  Mass. 

ATTN!  Dr.  Philip  Newman,  CRRK 

Mr.  N.  C.  Gerson 

Trapelo  Road 

South  Lincoln,  Mass. 

Dr.  Richard  B.  Barrar 
Systems  Development  Corp. 
2b00  Colorado  Avenue 
Santa  Monica,  California 

Columbia  University  Hudson  Laboratories 
P.O.  Box  239 

mS  Palisade  Street,  Dobbs  Ferry,  N.  Y. 
ATTN!  Dr.  N.  W.  Johnson 

Institute  of  Fluid  Dynamics 
and  Applied  Mathematics 
University  of  Maryland 
College  Park,  Maryland 
ATTN:  Dr.  Elliott  Montroll 

Department  of  Electrical  Engineering 

Washington  University 

Saint  Louis  $,   Mo. 

ATTN:  Professor  J.  Van  Bladel 

Department  of  the  Navy 

Office  of  Naval  Research  Branch  Office 

1030  E.  Green  Street 

Pasadena  1,  California 

Brandeis  University 
Waltham,  Mass. 
ATTN:  Library 

General  Electric  Company 
Mlctowave  Laboratory 
Electronics  Division 
Stanford  Industrial  Park 
Palo  Alto,  California 
ATTN:  Library 


Smyth  Research  Associates 
3555  Aero  Court 
San  Diego  3,  California 
ATTN!  Dr.  John  B.  Smyth 

Electrical  Engineering 
California  Institute  of  Technology 
Pasadena,  California 
ATTN:  Drfieorges  0.  Weill 

Naval  Research  Laboratory 

Washington  25,  D.  C. 

ATTN:  Henry  J.  Passerinl,  Code  5278A 


Dr.  George  Kear 
5  Culver  Court 
Orlnda,  California 

Brooklyn  Polytechnic 

85  Livingston  Street 

Brooklyn,  New  York 

ATTN!  Dr.  Nathan  Marcuvltz 

Department  of  Electrical  Engineering 
Brooklyn  Polytechnic 
85  Livingston  Street 
Brooklyn,  New  York 
ATTN:  Dr.  Jerry  Shmoys 

Department  of  Mathematics 
University  of  New  Mexico 
Albuquerque,  New  Mexico 
ATTN:  Dr.  I.  Kolodner 

Mathematics  Department 
Polytechnic  Institute  of  Brooklyn 
Johnson  and  Jay  Street 
Brooklyn,  New  York 
ATTN!  Dr.  Harry  Hochstadt 

Ballistics  Research  Laboratory 
Aberdeen  Proving  Grounds 
Aberdeen,  Maryland 
ATTN:  Dr.  Pullen  Keats 

Dr.  Lester  Kraus 
1:935  Whitehaven  Way 
San  Diego,  California 

University  of  Minnesota 
Institute  of  Technology 
Minneapolis,  Minnesota 
Attn:  Dean  Athelston  Spilhaus 

Ohio  State  University 
Columbus,  Ohio 
Attn:  Prof.  C.  T.  Tal 
Department  of  Electrical  Eng. 

Naval  Research  Laboratories 

Washington  25,  D.  C. 

Attn:  W.  S.  Ament,  Code  5271 

Naval  Research  Laboratory 
Washington  25,  D.  C, 
ftttn:  Dr.  Leslie  G.  MoCracken,  Jr. 
Code  3933A 

Office  of  Naval  Research 
Department  of  the  Navy 
Attn:   Geophysics  Branch,  Code  hl6 
Washington  25,  D.  C. 

Office  of  Chief  Signal  Officer 
Signal  Plans  and  Operations  Division 
Attn:  SIGOL-2,  Room  20 
Com.  Liaison  Br.,  Radio  Prop,  Sect. 
The  Pentagon,  Washington  25,  D.  C. 

Defence  Research  Member 
Canadian  Joint  Staff 
2001  Connecticut  Street 
Washington,  D.  C. 

Central  Radio  Prop.  Lab. 
National  Bureau  of  Standards 
Attn:  Technical  Reports  Library 
Boulder,  Colorado 

D.  S.  Weather  Bureau 
U.  S.  Department  of  Commerce 
Washinpton  25,  D.  C. 
Attn:  Dr.  Harry  Wexler 
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Federal  Cormnunlcatlons  Commission 

Washin.'ton  25,   D.  G. 

Attn:  Mrs.  Barbara  C.  Crimea,  Librarian 

Upper  Atmosphere  Research  Section 
Central  Radio  Propagation  Laboratory 
National  Bureau)  of  Standards 
Boulder,  Colorado 

Argonne  Matlonal  Laboratory 

P.O.  Box  299 

Lemont,    Illinois 

Attn:     Dr.  Hoylande  D.  Young 

Bell  Telephone  Labs. 

Murray  Hill,   New  Jersey 

Attn:     Dr.  S.  0.  Rice,   3B  -  203 

Carnegie  Institute  of  Washington 
Dept,  of  Terrestrial  Magnetism 
52lil  Broad  Branch  Road,  N.  W. 
Washington  1?,  D.  C. 
Attn:  Library 

Georgia  Tech  Research  Institute 
225  N.  Avenue,  N.  W. 
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